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Abstract 

In a previous paper of ours [13| we have considered the incompressible 
Navier-Stokes (NS) equations on a d-dimensional torus T , in the functional 
setting of the Sobolev spaces ]HI"(,(T°') of divergence free, zero mean vector 
fields (n > d/2 -f- 1). In the cited work we have presented a general setting 
for the a posteriori analysis of approximate solutions of the NS Cauchy prob- 
<^ I lem; given any approximate solution Ua, this allows to infer a lower bound Tc 

C\J ■ on the time of existence of the exact solution u and to construct a function 

^ . IR„ such that \\u{t) — Ua(t)||n ^ ^n(i) for all t G [0,Tc). In certain cases 

pg \ it is Tc = +c«, so global existence is granted for u. In the present paper 

,^ ■ the framework of [H] is applied using as an approximate solution an expan- 

O . sion u^{t) = '^j=QR-'uj{t), where R is the Reynolds number. This allows, 

C^ I amongst else, to derive the global existence of u when R is below some crit- 

ical value i?* (increasing with A^ in the examples that we analyze). After a 
general discussion about the Reynolds expansion and its a posteriori analysis, 
^ I we consider the expansions of orders A = 1, 2, 5 in dimension d = 3, with the 

initial datum of Behr, Necas and Wu |lj. Computations of order A = 5 yield 
a quantitative improvement of the results previously obtained for this initial 
datum in [Ijj, where a Galerkin approximate solution was employed in place 
of the Reynolds expansion. 

Keywords: Navier-Stokes equations, existence and regularity theory, theoretical 

approximation. 
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1 Introduction 

The incompressible Navier-Stokes (NS) equations with no external forces and peri- 
odic boundary conditions can be written as 

■£ = uAu + ?{u,u) , (1.1) 

where: z/ G (0, +oo) is the viscosity coefficient; u = u{x, t) is the divergence free 
velocity field; the space variables x = {xs)s=i,...,d belong to the torus T*^ (and yield 
the derivatives ds := d/dxg)] A := J2s=i ^ss is the Laplacian. Furthermore, 7 is the 
bilinear map defined as follows: for all sufficiently regular velocity fields v,w on T'^, 

T{v,w) := -£{v.dw) (1.2) 

where {v»dw)r := ^^=1 VgdsWr (r = 1, ..., rf), and Sj is the Leray projection onto the 
space of divergence free vector fields. The dimension d is arbitrary in the general 
setting of this paper, but we put (i = 3 in the application of the last section. From 
the physical viewpoint, the velocity fields (at any fixed time) should be maps from 
T"^ to R*^; however, we can harmlessly consider maps T'^ — )■ C^. 

Let us introduce the rescaled time t and the Reynolds number i?, setting 

t := z/t , R= - ■ (1.3) 

then Eq (II. ip takes the form 

^ = AM + i?T(M,M) , (1.4) 

ot 

that will be the standard of this paper. Our functional setting for Eq. (I1.4p is based 
on the Sobolev spaces 

Ml,{T^) = ni, ■= {v.T'^^ C^ \{v) = 0,divw = 0, V^A''\ e L'(T'^)} (1.5) 

(where ( ) indicates the mean over T*^); for any real n, the above space is equipped 
with the inner product {v\w)n '■= (y— A v\y/—A w)l2 and with the corresponding 
norm || ||„. 

In our paper [14J we have outlined a general framework to obtain quantitative 
information on the exact solution u of the NS Cauchy problem analyzing a posteriori 
an approximate solution. To be more precise, consider the NS equation (II. 4p with 
a specified initial condition u{x,0) = m*(x); let Ma : T"^ x [0,Ta) — )■ C^ be an 
approximate solution of this Cauchy problem, and consider (for n > d/2 + 1) the 
Sobolev norms 



dzi 

\\{-g^ - Au^- R?{u^,U^)){t)\\n , \M0)-U4n, (1-6) 

ll^aWIU, ||na(t)||„+l, (1.7) 

where t G [0,Ta) and u^{t) := Us,{-,t). The norms in fll.6p control the differential 
and datum errors of Ma, while the norms in (11.71) refer to the growth of u^- The 
approach of [11] relies on the so-called control inequalities; these consist of a differ- 
ential inequality and of an inequality on the initial value, determined by the norms 
(I1.6P (I1.7P and involving an unknown function 3^„ : [0,Tc) — )■ [0, -|-oo). Assume 
the control inequalities to have a solution !3^„, with a suitable domain [0,Tc); then, 
according to [14|, the solution u of the NS equation (ll.4p with initial datum u^, 
exists (in a classical sense) on the time interval [0,Tc), and its distance from the 
approximate solution admits the bound 

\\u{t) - U,{t)\\n ^ %n{t) for t e [0,re) . (1.8) 

(For similar statements on the NS equations or other nonlinear evolutionary PDEs, 
see [3J [lOj LH] [12] [H].) 

In the present paper we apply the above framework choosing as an approximate 
solution a polynomial in R of the form 

N 

u\t):=Y,R'u,{t), (1.9) 

j=0 

where the terms Uj{t) are determined requiring the differential error to be 0{R^^^) 
for i? — )• 0. We emphasize that, in our approach, the order N could be large but 
is fixed; so we are not considering the A^ — )■ -|-oo limit, i.e., the solution of the 
Cauchy problem for Eq. (11.40 via a power series u{t) = X]j=o R^Uj{t). A theoretical 
analysis of the convergence issue for such a series, in suitable function spaces, has 
been developed by some authors, especially Cannone P] and Sinai [TH]. However, 
the approaches of these authors yield convergence conditions (local or global in 
time) which depend on the initial datum u^ only through its norm; moreover, these 
authors have dedicated little attention to the strictly quantitative aspects of their 
analysis (such as the evaluations of the constants in certain inequalities). On the 
contrary, our approach based on a finite order approximant u^ as in (II. 9p has the 
following features. 

(i) We produce estimates on the interval of existence of the exact solution u (and 
on its distance from u^) which depend on the fine structure of the initial 
datum M^, and not only on its norm; the specific features of the initial datum 
yielding these estimates are encoded in the expression of the differential error 
du^/dt - Au^ - RT{u^,u^). 



(ii) Our analysis is fully quantitative: it relies on explicit expressions for u^ and 
its errors and uses, amongst else, the estimates of [13] [15] on the constants in 
certain inequalities about 7. 

The above setting invites a computer assisted approach: this can be readily set 
up when the NS initial datum m* is sufficiently simple, say, a Fourier polynomial. 
In this case the approximant u^ and the errors du^/dt — Au^ — i?T(-u^,M^), 
u^{0)—u^:, with their Sobolev norms, can be determined via any package for symbolic 
computation; after this a solution 3?„ for the control inequalities can be obtained 
numerically. More precisely, one can try to satisfy them as equalities: this amounts 
to solve the Cauchy problem for a simple ODE in the unknown function t i— )■ !3i„(t). 
This "control Cauchy problem" is easily treated numerically. 

In the present paper the above procedure is described in general terms and then 
applied with d = 3 and n = 3, choosing for u^, the so-called Behr-Necas-Wu (BNW) 
initial datum [Ij. For the practical implementation we use MATHEMATICA on a 
PC; first we work with A^ = 1 to introduce the method, and then pass to the orders 
N = 2, N = 5. 

In all the above cases, the control Cauchy problem has a solution 3^3 of domain 
[0, +oo) if R is below some critical value i?* (depending on A^); in this situation we 
can grant global existence for the solution u of the NS Cauchy problem, and the 
inequality \\u(t) — 'U^(t)||3 ^ '^sit) holds for all t G [0, +oo). For R above i?*, the 
control problem has a solution Jl^ on a bounded interval [0,Tc); so, the existence of 
the solution u of the NS Cauchy problem and the inequality \\u{t) —u^(t)\\3 ^ 3^3 (t) 
are granted only on this interval. Passing from A^ = 1 to A^ = 2, and from N = 2 
to N = 5, the critical value R^, increases. 

In the final part of the paper the outcomes of the above Reynolds expansions are 
compared with the results of [2] , where an approximate NS solution was constructed 
for the BNW initial datum using the Galerkin method with a set of 150 Fourier 
modes; it turns out that this Galerkin approach is quantitatively equivalent to the 
Reynolds expansion of order N = 2, while the expansion of order A^ = 5 gives 
much better results concerning the global existence of u and its distance from the 
approximate solution. 

The paper is organized as follows. After fixing some basic notations. Section [2] 
reviews the general setting of [2] for approximate NS solutions (in a reformulation 
suitable for our purposes, where the Reynolds number R = l/u and the rescaled 
time t = ut are preferred to the variables u, t). Sections |3] and Hj present the general 
Reynolds expansion (II. 9p . the control inequalities for it and some basic computa- 
tional rules for the practical implementation of this approach. Section \5\ applies the 
previous framework to the the BNW initial datum, and makes a comparison with 
the Galerkin approach of [H] for the same datum. 

In a forthcoming paper [16] the Reynolds expansion and the related control equa- 
tion will be applied to the BNW datum for larger values of A^, and then employed 



for other initial data of interest in the NS community, namely, the "vortices" of 
Taylor-Green [20] and Kida [6] . 

2 Preliminaries 

Throughout the paper we fix a space dimension d G {2,3,...}; in the application 
of section [5] we will put d = 3. For a, h in R*^ or C^ we put a»h := Ylir=i ^rK and 
\a\ := \/a»a, where ~ denotes the complex conjugation. 

Let us consider the li-dimensional torus T"' := (R/27rZ)'^; the setting outlined 
hereafter for function spaces and NS equations on the torus is similar, but not 
identical to the one of [12] [H]. A main difference between these papers and the 
present one is that here we deal with spaces of complex valued, rather than real 
valued functions; some results in [12] [13] (and in other works of ours), originally 
formulated in a real framework, hold as well in the present, complexified setting. 

In the sequel we employ the Fourier basis made of the functions (p) 

Cfc : T^ ^ C , efc(x) := e'''"' {k G Z^) . (2.1) 

Let us consider the space D'iT'^, C) = D' of complex distributions on T*^, i.e., the 
topological dual oi D := C°°(T'^, C) (p). Any v G D' has a weakly convergent Fourier 
expansion v = ^^^i^dVuek, where Vk '■= (27r)~'^(t>, e_fc) (the notation indicates the 
action of v on the function e_fc G D). 

We will often be interested in the spaces L^(T'^, C) = L^ and, in particu- 
lar, in the complex Hilbert space L^ with the product {v\w)l'2 := Jrj,ddxvw = 
(27r)'^X]feezd^'^fc and the corresponding norm || 11^2. For all n G R, the n-th 
Sobolev space of zero mean functions on T*^ is 

if^(T'^) = H^:= IveD' \ {v) = , V^A^v G L^ } (2.2) 

= LeD'\vo = 0, Yl I'^r^l^fcT < +00 } 

k£Z<i\{0} 

(in the above (v) G C is the mean of v, i.e., by definition, the action of v on the test 
function l/(27r)'^; moreover, \/— A v := J^kez^uo} m^'^k^k)- H^ is a Hilbert space 
with the inner product and the norm 

{v\w)^ := {V^\\V^K"w)l^ = {27TyJ2 l^r^^^fc. II^IU := ^H^ ; (2-3) 

kez^xio} 



^In [12] [M], the normalization factor (27r)^''/^ was included in the definition of Cfc. 

•^The space of complex distributions carries a complex conjugation , see, e.g., jl2] ; real distri- 
butions V on T'' are characterized by the condition v = v. To avoid confusion, we mention that 
in [12] [14] the space of complex distributions on T'' is indicated with D^, while the symbol D' 
is employed for the subspace of real distributions. Similar comments could be done about the 
Sobolev spaces mentioned later. 



iim^n, then HJ^ C H^. 



The Laplacian and its semigroup. Let us consider the operator A := X]rs=i ^^s • 
D' — )■ D'; of course Ae^ = —\k\'^ek for all k G Z"^. For each n G R, A maps 
continuously if""*"^ into H^, with ||At>||„ ^ 1 1 "^1171+2 for all v G if^"*"^. We can define 
a semigroup (e*^)jg[o,+oo) of linear operators on D', putting 



e'^:D'^D', v^ e'\ := J^ e~''''*^fcefc . (2.4) 

Let n G R. The following holds: 

e'^H^ C H;^ , \\e'\\\n < e-'\\v\\n for t G [0, +oo), t; G if ^ ; (2.5) 

gtA^n-i ^ ^n . g^ ^ ^ij^^^g^ ^^^j^ j^^, (independent of n) such that (2.6) 
||e*^t;|U ^ ;u(t)||t;||„_i for t G (0, +oo), t; G //„"-' 

(for the proof of (12.61) see, e.g., [12], that also gives an explicit expression for fi 
implying /j,(t) = 0(1/ \/t) for t — )■ 0^ and /i(t) = e~* for t large). The map {t,v) i— )■ 
e*^f is continuous from [0, +oo) x H^ to if" and from (0, +oo) x H^~^ to H^. 
Moreover, 

^ G H:^+^ ^ (t ^ e*^t;) G C([0, +oo), H:;+^) n C^([0, +oo), H^) , (2.7) 

^(e*^t;) = A{e'\) for t G [0, +oo). 

To go on, let 

/GC([0,+oo),ii,^i); (2.8) 

for each t G [0, +cx)) the function s G (0,t) h> e(*-")^/(s) is in L\{0,t),H;^+^) 
because, on the grounds of (12.61) . ||e'^*~'^^^/(s)||„+2 ^ /i(t — s)||/(s)||.„+i; therefore, 
the definition 

F{t) := f dse^'-'^^fis) for t G [0, +oo) (2.9) 

Jo 

produces a function F G C([0, +oo), if""^^). This function is also in C*^([0, +cx)), H"), 

with 

dF 

— it) = AF{t) + f(t) for t G [0, +oo) . (2.10) 

Vector fields on T'^. Here and in the sequel, "a vector field on T'^" means "a 
C^- valued distribution on T'^". We write ©'(T'^) = D' for the space of such distri- 
butions; these can be identified with d-tuples v = (t>i, ..., f^), where Vr G D' for each 
r. Partial derivatives, the Laplacian A and the operators \/—A , e*^ are defined 
componentwise as maps from D' to D'. Any f G ©' has a weakly convergent Fourier 
expansion v = J^keZ'^'^k^k, with coefficients Vk G C^. 



In the sequel IJ'{T'^) = LP denotes the space of L^ vector fields T'' — )■ C^. For 
each n G R, the n-th Sobolev space of zero mean vector fields on T'^ is 

H;^(T'^) = W^ := IveB' \ (v) = 0, y^A\ G V } (2.11) 

= Iv eB' \vo = , Y^ |A;p"|t;fc|2 < +00 } 
/cez<*\{o} 

(in the above, the mean {v) G C^ is defined componentwise; we note that H" = 
{v = {vi, ...,Vd) G D' I f r G if" for each r}). H" is a Hilbert space with the inner 
product and the norm 

{v\w)n:={V^A^v\V^K"w)L^ = {2nyJ2\^\''''^''^k, ||w|U := V(^|^)n • (2.12) 

A:GZ<*\{0} 

Eqs. (I2.5p - (l2.10p and the related statements have obvious analogues, where if" is 
replaced by H" for any n. 

Divergence free vector fields; the Leray projection. The space of divergence 
free vector fields on T'^ is 

B'^ := {v eB' \ divw = 0} = {v eB' \ k.Vk = V/c G Z'^ } . (2.13) 

The Leray projection is the linear, surjective map 

£ : D' ^ D'j, , v^ 2v:=Y^ {S.kVk)ek ; (2.14) 

fcGZ'* 

here £,k is the orthogonal projection of C^ onto the orthogonal complement of k, 
i.e., 

£oC = c, S.kC = c--^k for c G C^, A; G Z"^ \ {0} . (2.15) 

m 

For each n G R, the n-th Sobolev space of zero mean, divergence free vector fields is 

H^„:=D',nH^ (2.16) 

this is a closed subspace of H", and thus becomes a Hilbert space with the restriction 
of the inner product ( | )«• One has 

2.W^=ni^, ||£t;||„ ^ ||t;||„ for n G R, t; G H^^ . (2.17) 

The spaces (I2.16P are the basis of our treatment of the NS equations; again, we have 
analogues of Eqs. (l2.5p -f l2rT0|) and of the related statements, where if" is replaced 
by HI"q for any n. 



The NS bilinear map. Consider two vector fields f , w on T'^ such that f G L^ 
and drW G J-? for r = 1, ...,(i; then we have a well defined vector field v»dw G L-*^ 
of components {y»dw)r := X]s=i "^st^s^^ (which has mean zero if divt> = 0); we can 
apply to this the Leray projection 2, and form the (divergence free) vector field 

y{v,w):=-2{v,dw) . (2.18) 

The bilinear map T: (f , w) i— )■ y(w, w), which is a main character of the incompress- 
ible NS equations, is known to possess the following properties. 

(i) For each n > d/2, 7 is continuous from M^^ x M'^'^^ to H^^,; so, there is a 
constant Knd = Kn such that 

mv,w)\\n ^ K^\\v\\n\\w\\n+i for V G H^„, w G e^+1 . (2.19) 

(ii) For each n > d/2 + 1, there is a constant Gnd = Gn such that 

\{?iv,w)\w)r,\ ^ ajvUHll for V G HS„, w G H^+i (2.20) 

(this result is due to Kato, see ^). 

From here to the end of the paper, Kn and Gn are constants fulfilling the previous 
inequalities (and not necessarily sharp). From [13] [15] we know that we can take 

Ks = 0.323 , Gs = 0.438 if t/ = 3 ; (2.21) 

these values will be useful in the sequel. 

The NS Cauchy problem. From here to the end of the paper, we fix a Sobolev 
order 

ne (- + l,+oo) . (2.22) 

Let us choose a Reynolds number Re [0, +oo) and an initial datum 

u, G mi^^ . (2.23) 



2.1 Definition. The (incompressible) NS Cauchy problem with Reynolds num- 
ber R and initial datum m^, is the following: 

Find u G C([0, T), e^^+2) p ^iQo, T), H^^J such that (2.24) 

du 

-— = Au + RT(u,u) , u(0) = u^ 

dt 

(with T G (0, +oo], depending on u). 



It is known [5J that the above Cauchy problem has a unique maximal (i.e., not 
extendable) solution; any solution is a restriction of the maximal one. (u) 

Approximate solutions of the NS Cauchy problem. Let us briefly rephrase 
the results of [H], with the notations of the present paper (this is advisable, since 
in [13] we wrote the NS equations in the form (II. ip . rather than (II. 4p ). We consider 
again the Cauchy problem (I2.24p . for given n, i?, u^ as above; the definitions and the 
theorem that follow are reported from [2], with obvious adaptations. 

2.2 Definition. An approximate solution of the problem Ii2.24\) is any map 
Ma e C([0, Ta), Ml+^) n Ci([0, Ta), H^J (with Ta G (0, +oo];. Given such a function, 
we stipulate (i) (ii). 

(i) The differential error of u^ is 

du 

-^-Au,-R7{u,,u^) GC([0,Ta),H^J ; (2.25) 

the datum error is 

Ma(0) -u,e B.l+^ . (2.26) 

(ii) Let 771 G R, ?Ti ^ 72. A differential error estimator of order m for u^ is a function 
em G C([0, Ta), [0, +oo)) such that 

du 
\\{^-Au,-R7{u,,u,)){t)\\m^emit) fortG[0,Ta). (2.27) 

Let m ^H, m, ^ n + 2. A datum error estimator of order m for Ma is a real number 
6m G [0, +oo) such that 

\\UM - U4rn ^ 5m ; (2.28) 

a growth estimator of order m for Ua_ is a function Vm G C([0,Ta), [0,+oo)) such 
that 

||«a(t)|U ^ ^m(t) forte [0,Ta) . (2.29) 

In particular the function emit) '■= \\{dus^/dt — Aua — RV^Ua, Ua)){t)\\m, the number 
6m '■= ll^a(O) — M*||m 3.nd the function Vmif) := ||Ma(^)||m "^111 be called the tautolog- 
ical estimators of order m for the differential error, the datum error and the growth 
of Us,. 

From here to the end of the section we consider an approximate solution u^ of the 
problem (I2.24p . with domain [0,Ta); this is assumed to possess differential, datum 
error and growth estimators of orders n oi n + 1, indicated with e„, 5„, P„, 'Dn+i- 



^It is known as well that, if u solves (|2.24p . the function {x,t) i— > u{x,t) is smooth on T'^ x 
(0,T) (see, e.g.. Theorem 15.2 (A) of 0); one could give stronger regularity results with suitable 
assumptions on u*. In the sequel we will not be interested in such regularity matters. 



2.3 Definition. Let !3^„ G C([0,Tc), [0, +00)), with T^ G (0,Ta]. This function is 
said to fulhll the control inequalities if 

d+'Jl 

_J1 ^ _3?^ + R{GnVn + Kr,Vn+i)X. + RGnX^ + £„ in [O, Tj, (2.30) 

3?„,(0) ^ 5„ . (2.31) 

In the above d'^/dt indicates the right, upper Dini derivative: so, for all t G [0,Tc), 
{d+'Jlnldt){t) := limsup;,_,o+[3^„(t + h) - 3^„(t)]//i. 

2.4 Proposition. Assume there is a function "Rn £ C{\^,Tc)., [0, +c>o)) fulfilling 
the control inequalities; consider the maximal solution u of the NS Cauchy problem 
l\2.24\) . and denote its domain with [0,T). Then 

T^Tc , (2.32) 

Mt) - u,{t)\\n ^ Oln{t) forte [0,Te) . (2.33) 

For completeness we mention that any function 3l„ G C([0,Tc),R) fulfilling the 
inequalities (12.301) (12.311) is automatically nonnegative (p). 

Paper [14J presents some applications of the previous proposition, where Ua is 
constructed by the Galerkin method. (For completeness we mention that the frame- 
work of [Ilj also covers the case of the Euler equations, i.e., the limit case z/ — )■ 
of (II. ip which is formally equivalent to i? — ?■ +00; some applications to the Euler 
equations have been considered both in [T^ and in |9].) 

In the sequel of this paper we present an application of Proposition 12. 4[ choosing 
for Ma a polynomial in R (see Eq. (ll.9p ). In the next two sections we develop this 
approach in general terms, giving the error estimators for an approximate solution 
of this kind; in the final section we apply this procedure choosing for m* the BNW 
initial datum. 

3 Reynolds number expansions as approximate 
NS solutions 

Let us recall that n G {d/2 + 1, +00), and consider the NS Cauchy problem (I2.24p 
with R G [0, +00) and datum m^, G IH["+^. Let us choose an order A^ G {0, 1, 2, ...} 



^In fact the zero function fulfills relations analogous to ()2.30|) ()2.31|) . with > replaced by ^; 
by standard comparison results, this implies 3?„(i) ^ for all t e [0,Tc). The comparison results 
required to prove the last statement have been presented in [14] (see Appendix A of the cited paper 
and references therein). 



and consider as an approximate solution for f l2.24p a polynomial of degree N in R, 
of the form 

TV 

u^ : [0, +00) ^ mi+^ , t ^ u^{t) := ^Wuj{t) , (3.1) 

j=0 

Uj e c([o, +oo)e;!+2) p ^iQo, +00), e^j for j = o, ..., n- 

the functions Mj herein are to be determined. 
3.1 Proposition, (i) Let u^ he as in ^3.1\) . The datum and differential errors of 

N 

u are 

N 

M^(0) -u, = (uo(0) -u,) + Y^ R^u,{0) ; (3.2) 

dt 

, AT , J j-1 2Ar+l TV 

= (^-Ano)+5^i?^[^-A«,-5^y(«,,«,_,-i)]- 5^ i?^5^ y(n,,n,_,_i). 

j = l £=0 j=N+l e=j~N-l 

(a) One can define recursively a family of functions Uj G C([0, +00), H^^^) fl 
C^{[0,+oo),Ml^) setting 

Uo{t) := e*^M* for t G [0, +cx)) , (3.4) 

^j(^)-=y, cise(*-^)^y(n,(s),M,_,_i(s)) /ortG[0,+cx)), j = l,...,iV. (3.5) 

lyzi/i i/izs choice of uq, ...,un the coefficients of R^,R^, ...,R'^ in Eqs. (E^ and 
^3.3\) vanish, so that 

M^(0)-M* = 0; (3.6) 

j„,7V 2JV+1 TV 

^ _ An^ - T(«^, n^) = - ^ i?^- ^ T(n,, n,_,_i) . (3.7) 

j=Af+l £=j-Af-l 

The second equation implies 

||(^-An^-i?y(n^«^))(t)|U (3.8) 

2N+1 N 

^KnY^ R^ Y^ \\ui{t)\\n\\uj.t-i{t)\\n+i for all t E [O, +oo) . 

j=N+i e=j^N+i 
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Proof, (i) Eq. ( 13 ■2p is obvious. Let us prove Eq. ( I3.3p : to this purpose, we note 
that 

A,,N H ^ N N 

(JjV (JjL 

3=0 i=0 h=0 

AT , N 

j=o e,h=o 

, AT , 2Af+l 

/at, := {(£,/i) e {0,...,iV}2 I £ + /i + l = j} . 
One easily checks that 

je{l,...,N} => /^^. = {(£,j--£_i) |£e{0,...,j-l}}, 

jG{iV + l,...,2iV + l} ^ /^^. = {(£,^-_£_l) \ie{j-N-l,...,N}}; 

this readily yields the thesis (13.31) . 

(ii) First of all (by an obvious vector analogue of (12. 7p ). Eq. (13.41) defines a function 

Mo e C([0,+oo),e^+2) nCi([0,+cx)),e2j such that 

uo(0) = u^ , — = Auo ; (3.9) 

to go on, let us prove by recurrence over j that the functions Uj are well defined and 
belong to C([0, +oo), H^+2) n C\[0, +oo), H^J, for all j ^ A^. For j = 0, the thesis 
has been already estabhshed. Now, let j G {1, ..., A^} and assume the thesis to hold 
up to the order j — 1; then the functions t G [0, +cxd) t-)- 'J'{ui{t),Uj_i_i{t)) {i = 
0, ...., j — 1) are in C([0, +oo),E["+-'^) due to the properties of CP. By the obvious vector 
analogues of the considerations accompanying Eqs. (12. 9 p (I2.10p . with / replaced by 
7{u£,Uj-i-i), we see that the functions t t-)- Jg (ise*^*~*^^T('U^(s),Mj_^_i(s)) belong 
to C([0, +oo)EI^+^) nC^([0, +oo), E["o); the same can be said of Uj, which is the sum 
over £ of these functions. 

Let us prove that, with the present choice of the functions Uj {j = 0, ..., A^), the 
coefficients of i?",...,/?^ in Eqs. (13.21) (13.31) are zero. First of all, the definitions 
(13. 4p (13. 5p imply Mo(0) = u* and Uj{0) = for j = l,...,iV; so, the coefficients of 
R^, ..., i?^ in (13. 2p are zero. To go on, we note that Eq. (13. 9p indicates the vanishing 
of the coefficient of i?° in (13. 3p . For j = 1,...,N, using Eq. (I2.10p (in a vector 
analogue, with / replaced by ^^CP(m^,Mj_£_i)), we see that 

^(t) = Au.it) + J]T(«,(t),«,_,_i(t)) ; (3.10) 

£=0 
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this proves the vanishing of the coefficient of W in Eq. fl3.3p . 

Now Eqs. (13. 6p (13. 7p are obvious. Eq. ([321) and the inequahty ||CP('U^(t),'Uj_£_i(t))| 
^ -ft'„||'U^(t)||„||'Uj-^-i(t)||n+i gives immediately Eq. (13. 8p . D 

Let A^ G {0,1,2,...}; define u^ and Uj {j = 0,...,N) via Eqs. (^ ([32D fl33D . 
We regard u^ as an approximate sohition of the Cauchy problem (12.240 . of domain 
[0, +oo); this has the tautological datum error and growth estimators 



5n := 0, V„{t) := ||n^'(t)||„, P„+i(t) := h''(t)|U+i ■ (3.11) 

The differential error du^ /t — Au'^ — y(u^, u^) can be expressed via Eq. ( 13.70 : we 
can use for it the tautological estimator 

enit) := II (^ - A«^ - i?T(w^, «^)) (t) lU (3.12) 

or the rougher estimator indicated by (13. Sp . i.e., the function 

2N+1 N 

enit) := Kn Y, ^' H h£(t)||nh,-£-i(t)|U+i . (3.13) 

j=N+l £=j-N+l 

Now, using Proposition 12.41 with Wa = u^ and the above estimators we obtain the 
following. 

3.2 Corollary. Let N e {0,1,2,...}; deGne 5n,Vn,Vn+i via Eq. (EUP and 
e„ via Eq. i3.12\] or / I3.J3)) . for t G [0, +oo). Suppose there is a function Jin € 
C([0, Tc), [0, +CX))), witli Tc G (0, +cx)], fulfilling tlie control inequalities ^TM) H^) . 
Consider the maximal solution u of the NS Cauchy problem A2.24\) . of domain [0, T); 
then 

T^T,, ||M(t)-M^(t)|U^3^„(t) forte[0,T,). (3.14) 

In particular, the control inequalities h2. 30j) / IZ31J) are satisfied by any function 
X„ G C\[0, Tc), [0, +oo)) fulfilling the control Cauchy problem 

-^ = -Jln + R{GnVn + KnVn+l)Jln + RGnJll + enOn[0,T,), ^^^(O) = 0. (3.15) 

3.3 Remarks, (i) Even though rougher than the estimator (I3.12p . the estimator 
( I3.13P is interesting because its computation is less expensive; this is a relevant fact, 
especially in applications with a large A^. 

(ii) Let us consider the estimator ( I3.12p . writing via ( 13. 7p the error therein; alterna- 
tively, let us use the estimator (I3.13p . In both cases, it is natural to write 



enit) = R^^'init) for t G [0, +cx)) , (3.16) 
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where e{t) is a suitable function; this is a polynomial of degree A^ in i? in the case 
fl3.13p . and the square root of a polynomial of degree 2N in R in the case (13.121) . 
Consequently, the solution of the Cauchy problem (13.151) can be written as 

3^„(t) = i?^+iX„(t) , (3.17) 

where Jin G C'^{[0,T^), [0, +oo)) is such that 

dJl 

" -Jin + RiGnVn + KnVn+l)Jin + R^'^'^GnJil + €„ , 3?„(0) = 0. (3.18) 



dt 



(iii) Eq. (I3.14p has a number of obvious implications; let us give two examples, 
(iiii) The inequalities ||u^(t)||„ - \\u{t) - u^(t)||„ ^ ||^(t)IU ^ Ik^Wlln + \\u{t) - 
u^{t)\\n, the definition of Vn{t) in fl3TIl) and Eq. (Km for Jin imply 

Vn{t)-Jin{t)^\\u{t)\\n^Vn{t)+Jin{t) for t G [0, T^) . (3.19) 

(iii2) For any f G M^, the equation ||f ||^ = (27r)'^^j^,g2d\{o| I^P'^l'^fcP implies (27r)'^/^|ffc| 
^ ||f ||n/|^|" for all k. This inequality with v = u{t) — u^ {t) and (I3.14p give 

(27r)^/>,(t) - uUt)\ ^ ^0- for keZ^\ {0} and t G [0, T,) . (3.20) 



4 Implementing the recursion relations (13.41) (13.51) 



The recursion relations mentioned in the title are the main characters of Proposition 
13. Ij their building blocks are, essentially, the linear maps W, /C defined as follows: 

W:i7^+2^C([0,+oo),i/,7+2)p^iQ0,+oo),i7,!^), z^Uz (4.1) 

where Uz : [0, +00) -^ HJ^^'^ , t ^ {Uz){t) := e'^z ; 

/C : C([0, +00), /f^i) ^ C([0, +00), if^^) n Ci([0, +00), if.7) , / ^ /C/ (4.2) 

where /C/ : [0, +00) -^ H^^^^ , t ^ (/C/)(t) := / ds e^'-'^^fis) 

Jo 

(as for the domains and codomains of the above maps, recall Eqs. (12.71) (I2.10p and 
the related comments). These maps have vector analogues, denoted for simplicity 
with the same letters, 

W:e;^+2^c([o,+oo),e[^+2)p^iQo,+oo),e^), z^Uz, (4.3) 
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these are defined rephrasing Eqs. (14. ip fl4.2p . In the sequel, it will always be clear 
from the context whether the symbols W, /C refer to the maps (14. ip (14. 2 p or to the 
maps gD (I12D. For z = (^1, ..., Zd) G H^f and / = (A, ..., /,) G C([0, +oo), H^^+i) 
we have 

Uz = {piz^,...Mzd) , /C/ = (/C/i,...,/C/,) . (4.5) 

With the above notations, for any u^, G IHI"j|^ the recursion relations (13.40 (13. 5 p can 
be written as 

Uo := Wm* , (4.6) 

i-i 
n, := Y, ^^(^^' «i-^-i) (j = 1' •••' N) ■ (4-7) 

Hereafter we give a set of elementary computational rules about the maps W, /C; 
these allow a straightforward computation of the functions Uo,Ui, ...,U]\f in Eqs. 
(14. 6 p (14. 7p . at least in the case when u^, is a Fourier polynomial (see the remark 
at the end of this section). Our computational rules involve the functions e^ G D, 
ek{x) := e^^'^ {k G Z'^) and 

5,,fe : [0, +CX)) ^ R , t H> 5,,b(t) := Te"'* (a, 6 G N) , (4.8) 

as well as the products 

Ba,b Ck : [0, +oo) ^D , t^ Te-'* e^ ; (4.9) 

they can be summarized in the following statement. 

4.1 Lemma. Let a, a', b, h' G N, fc, k' G Z*^ and s G {1, .., d}. Then 

ds{Ba,bek) = iksBa,bek ; (4.10) 

{Ba,bGk){Ba'^b' Gki) = Ba+a',b+b' ^k+k' ■ (^-H) 

In addition, let k ^ Q (so that the functions below are in the domains ofU or fC). 
Then 

Uck = 5o,|fc|2efc ; (4.12) 

K(fi„,ltp et) = ^°+'-"''' et . (4.14) 
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Proof. Eqs. (14.1 Op fl4.1ip are obvious consequences of the definitions of e^ and Ba^- 
Eq. fl4.12p is just a reformulation of the relation e^^Ck = e"'*^' * Cfc. 

In order to derive Eqs. (14.130 fl4.14p . we note that (both for b ^ |fcp and for 

b = i^n, 

/C(S,,feefc)(t) = / dse^'-'^'^is'^e-'^ek) (4.15) 

Jo 



Let h 7^ |A;p. Then, a change of variables a = {h — \k\'^)s gives Jg ds s°'e ^^ ''^' •** 
= 1/(6- |A;|2)"+i/p(''-l'l')*rfaaV-- thus 



e I'^P* 



/C(5„,, efe)(t) = ^^_^j^^,^,+, 7(« + !,(&- I A:nt) e, , (4.16) 

where we have introduced the incomplete Gamma function 

ry 
j{a,y):= daa'^-^e"" for y e R, a G {1, 2, 3, ...} . (4.17) 

Jo 

It is known that 

a £ 

j{a + l,y) = a\(^l-e-yJ2j^) (4-18) 

i=0 

(see, e.g., |T^, page 177, Eqs. (8.4.7) and (8.4.11)). Inserting this result into (I4.16P 
we obtain that, for b ^ |/i;p, 

(— Ifcpt '^ j.£ —bt \ 

(F^m=^-E(j3^pj^«rj«Je.; (4.19) 

this proves Eq. (I4.13p . Finally, in the case 6 = |/i;p we obtain from (I4.15P that 

/C(5„,|,pefc)(t) = e-l^l * U dss'^j e, = — -^ e^ ; (4.20) 

this proves Eq. fj4?T4|) . D 



4.2 Remark. Let us return to the recursion rules (14. 6 p (14. 7p . assuming that the 
initial datum u^ is a Fourier polynomial: by this we mean that u* has finitely 
many non zero Fourier coefficients. In this case, due to Eqs. (14. 6 p and (I4.12p all 
components of Mq are linear combinations of finitely many functions of the form 
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i?o,|fc|2 Cfc. Using Eq. fl4.7p with the resuhs of the previous lemma we see that, 
for j = 1,...,A^, each component of Uj is a finite hnear combination of functions 
of the form Ba,b^k- We note that each term CP('u^,-Uj_£_i) in Eq. (14. 7p can be 
calculated as follows: first of all, one computes each component of Ui»duj_£_i using 
elementary considerations of bilinearity, together with Eqs. fl4.10p f l4.1ip : next, one 
obtains CP(m^, Mj_£_i) = — £(%«9mj_£_i) using the expression ( l2.14p ( ETT^ for the 
Leray projection. 

5 An application with the Behr-Necas-Wu (BNW) 
initial datum 

Throughout this section we work with 

d = 3, n = 3 (5.1) 

and any Reynolds number R G [0, +oo). The Cauchy problem f l2.24p takes the form 

Find u e C([0,T),H^J n C^([0,T),H^J such that (5.2) 

du . , 

— = A'U + i?T(M,M) , u{fi)=u^. 

(JjL 

The initial datum m* in ]H[|g (in fact, in H^^ for any real m) is chosen as follows: 

U:tixi, X2, xs) := 2(cos(xi + X2) + cos(xi +X3), (5.3) 

— COs(xi + X2) + C0S(X2 + X3), — COs(xi + X3) — C0S(X2 + Xs)) . 

This is a Fourier polynomial; indeed 

u* = ^ u^k^k , (5.4) 

fc=ita,it6,±c 

a:=(l,l,0), 6:=(1,0,1), c := (0, 1, 1) ; 
M*,±a := (1,-1,0) , M*,±b := (1,0, -1) , M=^,±c := (0,1,-1) . 

This initial datum has been introduced by Behr, Necas and Wu in [1]; these authors 
have considered the datum 05.30 (15. 4p as the origin of a possible blow-up for the 
Euler equations (i.e., for the zero viscosity limit of (II. ip ). 

We have disputed the BNW blow-up conjecture in [8]; in the present paper, 
independently of any opinion on the validity of this conjecture, we consider the 
datum (15.31) (15. 4p in presence of viscosity and analyze it by the Reynolds expansion 
method of sections [3] and HI In the final subsection we compare the outcomes of 
this approach with the results on the BNW datum obtained for the viscous case in 
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], where a Galerkin approximate solution was employed in place of the Reynolds 
expansion. 

Setting up the Reynolds expansion. The expansion of order A^ relies on the 
function u'^ := J2j=o^'^'^j '■ [0, +00) — > ]HI|q where Uo,...,Uiy are computed via the 
recursion rules (14.61) (14. 7p . starting from the BNW datum (15. 4p : after finding u^ , 
one computes the related growth and error estimators ©3, 2)4,63 and sets up the 
control Cauchy problem 

_i = _3i3 + R^G^Vs + 7^3^4)3^3 + RGsOil + €3 on [0, T,), 3^3(0) = , (5.5) 
with K^ and G3 as in (I2.2ip . In the sequel 

3^3: [0,Tc)-> [0,+oo) (5.6) 



always denotes the maximal solution of this problem (jj). 

The above conceptual scheme has been implemented in the cases A^ = 1,2,5, 
on which we report in the sequel. The case iV = 1 is useful to illustrate in full 
detail the method; the other ones give more accurate results, at the price of more 
expensive computations. In all cases, for the practical computation of mo,mi, ••• and 
of the estimators 2)3,1)4, £3, we have employed MATHEMATICA in the symbolic 
mode. Then the control Cauchy problem has been solved numerically (using again 
MATHEMATICA) for several sample values of R; these numerical calculations are 
very realiable, since they concern a simple one- dimensional ODE. The outcomes of 
such computations give evidence for the following picture. 

(i) There is a critical Reynolds number /?* (depending on A^) such that T^ = +00 
for ^ -R ^ i?=K, and T^ < +00 for R > R^,. Moreover, for ^ R ^ R^, one has 

^3(t) -^ 0+ for t -^ +CX), while for R> R^ one has ^^it) — )■ +00 for t — )■ T^ . 

On the grounds of our general setting, the results (i) on the control Cauchy problem 
yield the following conclusion. 

(ii) Let us consider the maximal solution u of the NS Cauchy problem (15. 2p and 
the Reynolds expansion (for A^ = 1,2 or 5), with its critical number R^. For 
^ R ^ R^, u is global and \\u{t) - M^(t)||3 ^ ^^{t) for all t G [0, +cx)). For 
i? > -R=K, we can only grant that the domain of u contains the interval [0, Tc), 
and that \\u{t) -u^{t)\\3 ^ Jlait) for all t G [0,Tc). 



^More precisely: one considers the real valued C^ functions fulfilling ()5.5p . which are automat- 
ically nonnegative (see the comment immediately after Proposition 12 .4^ : the maximal C^ solution 
(i.e., the solution with the largest domain) is the function in (|5.6p . 
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In the sequel we give more detailed information about each one of the cases N = 
1, 2, 5. As expected, when A^ is increased the critical value R^, increases as well. 
Our tests also indicate that, for a given R, when N increases the same happens of 
Tc] on the contrary 3?3 becomes smaller, thus giving a more stringent estimate on 

\\uit)-u^it)h. 

Case AT = 1. This relies on 



u^ ■= uq + Rui : [0, +oo) -^ Ml 



5 

so 



(5.7) 



The expressions of Uj = (mj,i,Mj,2, ^j.s) for j = 0, 1 are as follows (with Ba^hif) '■'■ 
t"e~''*, as in the previous section): 

Wo,i = -Bo,2(e(_i -1,0) + e(-i,o,-i) + e(i,o,i) + e(i,i,o)) , (5.^ 

Mo,2 = -Bo,2(-e(-l -1,0) + 6(0,-1,-1) + 6(0,1,1) - 6(1,1,0) , 

■"0,3 = -Bo,2(-e(-i,o,-i) - 6(0,-1,-1) - 6(0,1,1) - 6(1,0,1)) ; 



2i 

UlA = Y5o,4(e(-2,-l,-l) + 6(-l,-2,-l) - 6(-i,-i,-2) + 6(1,1,2) " 6(1,2,1) - 6(2,1,1)) (5.9) 

2i 

+ -;r-D0,6(— 6(-2,-l,-l) — 6(-l,-2,-l) + 6(-i,-i,-2) — 6(1,1,2) + 6(1,2,1) + 6(2,1,1)) , 
^1,2 = ^-D0,4( — 6(_2,-l,-l) — 6(-i,-2,-l) — 6(-i,-i,-2) + 6(1,1,2) + 6(1,2,1) + 6(2,1,1)) 

2i^ . . 

+ V-"0.6l6(-2,-l,-l) + 6(-i,-2,-l) + 6(-i,-i,-2) - 6(1,1,2) " 6(1,2,1) " 6(2,1,1)) , 

2i 

Ul,3 = y5o,4(-6(_2,-l,-l) + e(-i,_2,-l) + e(_i,_i,_2) - 6(1,1,2) " 6(1,2,1) + 6(2,1,1)) 

2i^ , , 

+ -^-D0,6(6(_2,-l,-l) — 6(-i,-2,-l) — 6(-i,-i,_2) + 6(1,1,2) + 6(1,2,1) — 6(2,1,1)) • 

The next step is to compute the differential error du^/dt — Au^ — R7{u^, u^) (which 
can be expressed via (13.71) ) and the (time dependent) norms 



©3 := llu^lls, P4 := \W\\i, 63 : = 
These are as follows: 



du^ 



Au^-Ry{u^,u^)\\3. (5.10) 



Pa = 4v^(27r)3/2 [^0,4 + 18i?2(5o,8 - 25o,io + ^0,12) 
V^ = 8v^ (27r)3/2 [^0,4 + 54i?2(5o,8 - 25o,io + ^0,12) 



1/2 



1/2 



(5.11) 
(5.12) 
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£3 = 8 v^ (27r)3/2i?2 ^1065(So,i2 - 25o,i4 + ^o.ie) (5.13) 

-|l/2 

+3872i? (So,i6 — 4i?o,i8 + 6i?o,2o — 4_Bo,22 + -60,24) 

The above functions determine the control Cauchy problem (15. 5p . The numerical 
solution of this problem for many sample values of R yields a picture as in items 
(i)(ii), pageHTl with a critical Reynolds number 

i?, e (0.08,0.09) . (5.14) 

In Boxes la-Id we consider the case R = 0.08, giving information on the following 
functions of time: the quantity (27r)'^/2|u^(t)| for the wave vector k = (1, 1,0); the 
estimators V^ and £3; the solution 3?3 of the control Cauchy problem, which is global. 
In Boxes 2a-2d we consider the analogous functions in the case R = 0.09, where "Jl^ 
diverges at T^ = 2.153... (P). Each one of these boxes (and of the subsequent ones) 
contains the graph of the function under consideration, and its numerical values for 
some choices of t. 

One immediately notices that the functions in boxes of the types (a) and (b) 
(i.e., the norm of the Fourier component (1,1,0) and the estimator D3) are very 
similar in these examples (and in all the subsequent ones), even from the quantita- 
tive viewpoint. What really makes the difference among these examples (and the 
forthcoming ones) are the differential error estimator 63 and the solution 3^3 of the 
control Cauchy problem, considered in type (c) and (d) boxes. 

Case AT = 2. This relies on the function 

u^ ■= uq + Rui + R^U2 : [0, +cx)) -^ H^q . (5.15) 

Of course, Uq^Ui are as in Eqs. fl5.8p fl5.9p : U2 = (m2,1) W2,2?'W2,3) has a more com- 
plicated expression, and in order to save room we only report its first component. 
This is 

U2,l = g5o,2 (-e(_i -1,0) - 6(^1,0-1) - 6(1,0,1) - 6(1,1,0)) (5.16) 

+ --Bo,4 (e(o,_2,0) - 6(0,0-2) - 6(0,0,2) + 6(0,2,0)) 
+ — -Bo,6 (— 6(-3,-2,-l) — 6(_3,_i,_2) + 6(-2,-3,-l) — 6(_2,-l,-3) — 6(_i,_3,_2) " 6(_i,_2,-3) 

+46(_i,_i,o) + 46(_i,o,_i) — 8e(o,_2,o) + 86(0,0,-2) + §6(0,0,2) — §6(0,2,0) 

+46(1,0,1) + 4e(i,i,o) — 6(1,2,3) — 6(1,3,2) " 6(2,1,3) + 6(2,3,1) — 6(3,1,2) — 6(3,2,1)) 



^Here and in the sequel, an expression like r — a.bcde... means that a.bcde are the first digits 
of the MATHEMATICA output in the computation of r. 
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+ q-Bo,8 (e(-3,-2 -1) + e(-3 -1 -2) — e(_2 -3,-1) + e(-2,-i -3) + ^(-i -3,-2) + ^(-i -2 -3) 

— 2e(_i _i^o) — 2e(_i^o -1) + 3e(o -2,0) — 3e(o,o,-2) — 3e(o,o,2) + 3e(o,2,o) 
-2e(i,o,i) - 2e(i_i,o) + 6(1,2,3) + 6(1,3,2) + 6(2,1,3) - 6(2,3,1) + 6(3,1,2) + 6(3,2,1)) 

+ ^-Bo,14 (-6(-3,-2,-l) - 6(_3,_l,_2) + 6(_2,-3,-l) " 6(-2,-l,-3) " 6(-l,-3,-2) " 6(_l,_2,-3) 
—6(1,2,3) — 6(1^3,2) — 6(2,1,3) + 6(2,3,1) " 6(3^1,2) — 6(3,2,1)) • 

The next step involves the differential error dv? / dt — ^v? — R7{v?' , u^) and the (time 
dependent) norms 

du^ 
Vs := 11^13, ^4 := \\u%, £3 - 11^ - ^«' - ^y(w', w')ll3 • (5.17) 



For example, one has 

1296 So,4 + 288i?2 (_B + 84 Bq. - 164 Bo.io + ^IBq 12) 



^.^^(2.P 



+R'^ (16 5o,4 + 1056 5o,8 - 4544 5o,io + 16317 ^0,12 - 29496 ^0,14 + 17680 Bo. 
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+ 6174^0,20- 8232^0,22 + 1029 5o,28)] • (5.18) 

The expressions of P4 and €3 are not reported. The former has a complexity similar 
to that of T'3, the latter is lengthier; in fact, €3 has the form (27r)^/^(^ • j, Cj,f,i?-^i?o,fe)^^^, 
where Cj^b are rational coefficients and the sum involves 48 pairs {j,b), with j e 
{6,8,10}.' 

Let us pass to the control Cauchy problem (15. 5p : in the present case, the picture 
of items (i)(ii), page [17] is realized with a critical Reynolds number 

i?, G (0.12,0.13) . (5.19) 

Boxes 3a-3d are about the case R = 0.12, and give information on the functions 
already chosen for the previous tables; one of them is the solution 3^3 of the control 
Cauchy problem, which is global. Boxes 4a-4d are about the analogous functions in 
the case R = 0.13, in which 3^3 diverges at Tc = 2.604.... 

Case AT = 5. This relies on the function 

5 
u^ := Y^ R^Uj : [0, +00) -^ H^,, . (5.20) 

i=o 

The terms Wo,Mi,M2 are as before; the functions ^3,^4 and M5 have expressions of 
increasing complexity, that we cannot reproduce here. We only mention that each 
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one of the components u^^i, u^^2i ^5,3 is a linear combination with rational coefficients 
of 1924 terms of the form iBa^Ck, with a G {0, 1, 2}; the wave vectors k appearing 
at least in one of the three components of U5 are 174. 

The differential error du^/dt — Au^ — R7{u^, u^) and the norms 

dii 

V,:=\\u%, V,:=\\u%, 63 := ||— - An^ - i?T(n^M5)||3 (5.21) 

dt 

have very lengthy expressions. Each one of the functions (15.211) has the form 
(27r)^/^(V . ^ jjCj,a,6-R"'-Ba,fe)^^^, where the Cj^a,b are rational coefficients and the sum 
involves finitely many triples (j, a, b) of nonnegative integers; these triples are 204 
in the cases of P3 and 1^4, and 1734 in the case of €3. 

As in the other cases, the final step is the control problem (15. 5p : the picture 
outlined by items items (i)(ii), page [17] is now realized with a critical Reynolds 
number 

i?, G (0.23,0.24) . (5.22) 

In Boxes 5a-5d and 6a-6d we give some information about the cases R = 0.23 and 
R = 0.24, respectively. Boxes 7a- 7d are about the case R = 0.12; this choice is 
considered as well in the next subsection, where it is used for a comparison between 
the Reynolds expansion and the Galerkin approach. 

Comparison with the Galerkin approach. In ^14j we have considered the 
NS Cauchy problem with the BNW initial datum, using the Galerkin approximate 
solution u'^ that corresponds to a finite set of Fourier modes G C Z^ \ {0}. 
This approach relies on the (finite-dimensional) Galerkin subspace 

mg, := {v G D'j^o \vk = ioT k^G} (5.23) 

= {^.'^fcCfc I Vk G C^, k»Vk = for all k G G} 

k&G 

and on the projection 

qj^ : D'^0 ^ H^o . v= Y. ^'^^'^ ^ ^^""^ ■= H '"^''^ ■ (^-24) 

A:eZ<^\{0} fceG 

The Galerkin approximate solution of the Cauchy problem (15. 2p (with the BNW 
datum) corresponding to the set of modes G is the unique function m*^ such that 

dvP 

M^GC^([0,+oo),eg,) , —— = Au^ + R^^y{u^,u^), M^(0) = *P^M, (5.25) 

{JiL 

(thus, u'^ solves a finite-dimensional Cauchy problem; one can show the existence of 
a global solution, of domain [0, +C)o), using the fact that the L^ norm is a decreasing 
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function of time). Indeed, in [13] we considered an equivalent formulation of f l5.25p 
based on the viscosity z/ and on the unsealed time t, related to R and to the present 
time variable t via Eq. (11. 3p : in the sequel we will rephrase the results of [13] in 
terms of the variables R and t. 

For a given G, u^ can be computed numerically; more precisely, one solves 
numerically a system of ODEs for the Fourier coefficients 7^ in the expansion 



u"" 



it) = Y,lk{t)ek . (5.26) 



keG 



One can specialize to u^ the general framework for approximate NS solutions; in 
particular, one introduce the tautological estimators of orders 3 or 4 for the growth 
and the differential error of m*^, i.e., the functions 

V,{t):=\\u^{t)h, V,{t):=\\u''{t)\U. (5.27) 

Mt) ■■= II (^ - Au^ - i?T(«^, u^)) (t) II3 (te [0, +00)) 

(see [H] for an explicit expression of 63 in terms of the Fourier coefficients 7jt). The 
datum error m'^(O) — m* is zero if G contains the Fourier modes ±a, ±6, ±c involved 
in Eq. (15. 4p . Assuming this, one can analyze the Galerkin approximate solution in 
terms of a control Cauchy problem 

—l = -Jl^ + R(G3Vs + K^V^)Jls + RG3Jll+?son[0,T,), 3^3(0) = , (5.28) 

whose (maximal) solution 3?„ G C^([0,Tc), [0,+oo)) can be computed numerically. 
In [14] we have employed a set G of 150 modes, of the following form: 

G:=SU-S; -S := {-k \ k e S} ; (5.29) 

S := {(0, 0, 2), (0, 1,-3), (0, 1, 1), (0, 1, 3), (0, 2, 0), (0, 2, 2), (0, 3, -1), (0, 3, 1), (0, 3, 3), 

(1, -3, -2), (1, -3, 0), (1, -3, 2), (1, -2, -3), (1, -2, -1), (1, -2, 1), (1, -2, 3), 
(1,-1, -2), (1,-1,2), (1,0, -3), (1,0,1), (1,0,3), (1,1, -2), (1,1,0), (1,1,2), (1,2, -3), 
(1, 2, -1), (1, 2, 1), (1, 2, 3), (1, 3, -2), (1, 3, 0), (1, 3, 2), (2, -3, -3), (2, -3, -1), (2, -3, 1), 
(2, -3, 3), (2, -2, -2), (2, -2, 2), (2, -1, -3), (2, -1, -1), (2, -1, 1), (2, -1, 3), (2, 0, 0), 
(2, 0, 2), (2, 1, -3), (2, 1, -1), (2, 1, 1), (2, 1, 3), (2, 2, -2), (2, 2, 0), (2, 3, -3), (2, 3, -1), 
(2, 3, 1), (2, 3, 3), (3, -3, -2), (3, -3, 2), (3, -2, -3), (3, -2, -1), (3, -2, 1), (3, -2, 3), 
(3, -1, -2), (3, -1, 0), (3, -1, 2), (3, 0, -1), (3, 0, 1), (3, 0, 3), (3, 1,-2), (3, 1, 0), (3, 1, 2), 
(3, 2, -3), (3, 2, -1), (3, 2, 1), (3, 2, 3), (3, 3, -2), (3, 3, 0), (3, 3, 2)} . 

Let us summarize the results arising from this choice of G. 
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( I ) There are indications for the existence of a critical Reynolds number i?^ such 
that: Tc = +00 and 'Jl^ vanishes at +00 if ^ i? ^ i?^, while Tc < +00 and 
^3 diverges at %iiR> R^. One has R^ E (0.13, 0.14) g). 

(n) Let us consider the maximal solution u of the NS Cauchy problem (15.21) . For 
^ -R ^ -R^, (?) grants that u is global and \\u(t) — 'u'^(t)||3 ^ 3^3 (t) for all 
t G [0, +00). For R > R^ (i) only grants that the domain of u contains the 
interval [0,fc), and that \\u{t) -M^(t)||3 ^ %{t) for all t G [0,fc). 

As an example, Boxes 8a-8d report the main results about the Galerkin approach 
for R = 0.12. Both the general picture (T) (n) and the results for R = 0.12 indicate 
a substantial equivalence between the Galerkin approach with the set G of (I5.29P 
and the Reynolds expansion of order N = 2. In fact: 

(a) The expansion of order N = 2 yields a picture similar to (T) (n), with a critical 
Reynolds number /?* G (0.12,0.13) (see items (i)(ii) of page [T7l and Eq. (15.19^ 
for i?*). 

(b) The error estimators €3, €3 and the solutions Jl^, 3^3 of the control equations 
for the N = 2 Reynolds expansion and for the Galerkin approach with G in 
(I5.29P are have similar orders of magnitude in the case R = 0.12, here used 
as a test to make comparisons: see Boxes 3c-3d and 8c-8d for more detailed 
information. 

If we pass to the N = 5 expansion, we find a significant improvement with respect to 
the results of the above Galerkin approach. Let us recall that, for N = 5, the critical 
Reynolds number for global existence is in the interval (0.23,0.24); moreover, if we 
use again the case R = 0.12 for a comparison, we see that the error estimators 63 
and the function CR3 of the N = 5 expansion are much smaller than the homologous 
function 63, Ji^ of the Galerkin approach: see Boxes 7c- 7d and 8c-8d. In particular, 
the ratio Ji:i{t)/%{t) is of order 10"^: so, the bound \\u{t)-u^{t)\\s ^ Ji^it) is much 
more stringent than the bound \\u(t) — u'^(t)\\3 ^ 3^3 (^)- 



®The estimate coming directly from [14] is 1/8 < i?* < 1/7, whence 0.125 < i?* < 0.143; this 
estimate has been refined to 0.13 < i?* < 0.14 by a supplementary run of the MATHEMATICA 
program for [14j . 
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12 3 4 


Box la 


. iV = 1, fl = 0.08: the function 


7(t) : = 


2,r)3/2|„l^_^_gj(t)|. One has 7(0) = 


22.27..., 


7(0.5) = 8.189..., 7(1) = 3.012..., 


7(1.5) = 


-- 1.108..., 7(2) = 0.4076..., 7(4) = 


7.466... 


X 10~^, 7(8) = 2.504... X 10~'^, 


7(10) = 


4.587... X 10"** . 
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12 3 4 


Box lb 


. AT = 1. fl = 0.08: the function 


-Osit) ■■= 


||M^(t)||3. One has O3 (0) = 154.3..., 


I>3(0.5) 


= 56.94..., 133(1) = 20.90..., 


I'3(l-5) 


= 7.683..., -03(2) = 2.826..., 


1^3(4) = 


0.05176..., 03(8) = 1.736... X lO"'^, 


153(10) 


= 3.180... X 10"'' . 
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Box Ic. iV = 1, fl = 0.08: the function 


esit). One has £3(0) = 0, £3(0.14) = 2.266..., 


£3(1) = 0.04605..., £3(2) = 1.296... X 10~*, 


£3(4) = 8.108... X 10""-" . 
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Box Id. N ^ 1, R = 0.08: the function 
Di^(t). This appears to be globally defined, 
and vanishing at +00. One has 3^3(0) — 
0, 3^3(1) = 9.858..., 3^3(1.5) = 11.66..., 
313(2) = 10.23..., 3^3(4) = 2.283..., K3(8) = 
0.04547..., 313(10) = 6.162... x 10~^ . 
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Box 2a. JV = 1, fl = 0.09: the function 


7(t) := (27r)3/2|„l^ i,o)(*)l- One has 7(0) = 


22.27..., 7(0.5) = 8.188..., 7(1) = 3.011..., 


7(1.5) = 1.107..., 7(2) = 0.4075..., 7(4) = 


7.465... X 10"^ . 
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Box 2b. JV = 1, fi = 0.09: the function 


T>3(t) := ([^^(t)!^. One has I>3(0) = 154.3..., 


I>3(0.5) = 56.99..., 133(1) = 20.90..., 


©3(1.5) = 7.684..., 133(2) = 2.826..., 


I>3(4) = 0.05176... . 


3 




2.5 


A 


2 


\ 


1.5 


\ 


0.5 


V 






12 3 4 




Box 2c. JV = 1, fl = 0.09: the function 


£3(4). One has £3(0) = 0, £3(0.14) = 2.868..., 


£3(1) = 0.05829..., £3(2) = 1.640... X 10"*, 


£3(4) = 1.026... X 10"^ . 
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Box 2d. JV = 1, ij = 0.09: the function 


3J3(t). This diverges for t ^ Tc = 


2.153.... One has 3?3(0) = 0, 3?3(0.5) = 


6.624..., 3^3(1) = 23.20..., 3^3(1.5) = 46.60..., 


3?3(2) = 176.0... . 
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Box 3a. JV = 2, fl = 0.12; the function 


7(t) := (2^)3/2|„J^_^ ^^(t)\. One has 7(0) = 


22.27..., 7(0.5) = 8.184..., 7(1) = 3.009..., 


7(1.5) = 1.107..., 7(2) = 0.4072..., 7(4) = 


7.459... X 10"3_ ^(g) ^ 2.502... x 10"'', 


7(10) = 4.583... X 10^** . 
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Box 3b. JV = 2, fl = 0.12: the function 


T>3{t) := ||ti^(t)||3. One has C3 (0) = 154.3..., 


I>3(0.5) = 57.10..., Xi3(l) = 20.88..., 


I'3(l-5) = 7.672..., 153(2) = 2.821..., 


I>3(4) = 0.05168..., r>3(8) = 1.733... X 10^^, 


E'3(10) = 3.175... X 10"^ . 
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Box 3c. AT = 2, fl = 0.12: the function 


£3(4). One has £3(0) = 0, £3(0.20) = 0.2759..., 


£3(1) = 0.01007..., £3(2) = 1.688... X 10-*, 


£3(4) = 5.653... X 10"** . 
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Box 3d. JV = 2, i? = 0.12: the function 


3?3(t). This appears to be globally defined. 


and vanishing at +00. One has 3?3(0) — 0, 


«3(1) = 4.699..., 3i3(1.6) = 6.854..., 


^3(2) = 6.368..., 3i3(4) = 1.466..., ^3(8) = 


0.02914..., .■R3(10) = 3.949... X 10"^ . 
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Box 4a. JV = 2, fl = 0.13: the function 


7(t) := (27r)3/2t„l^ i,o)(*)l- 0"° ^'^'^ ^W = 


22.27..., 7(0.5) = 8.183..., 7(1) = 3.009..., 


7(1.5) = 1.106..., 7(2) = 0.4071..., 7(4) = 


7.457... X 10"^ . 
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Box 4b. JV = 2, fl = 0.13: the function 


I>3(t) := \\u^(t)\\3. One has I>3(0) = 154.3..., 


03(0.5) = 57.16..., Xi3(l) = 20.89..., 


I>3(1.5) = 7.671..., ■03(2) = 2.821..., 


r>3(4) = 0.05166... . 
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Box 4c. JV = 2, fl = 0.13: the function 


£3(4). One has £3(0) = 0, £3(0.20) = 0.3508..., 


£3(1) = 0.01280..., £3(2) = 2.146... X 10-*', 


£3(4) = 7.187... X 10-** . 
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Box 4d. JV = 2, fl = 0.13: the function 


K3(t). This diverges for t ^ T^ = 2.604... . 


One has ^3(0) = 0, 3^3(0.5) = 1.735..., 


^3(1) = 10.20..., 3?3(1.5) = 20.68..., 


^3(2) = 36.30..., 3?3(2.5) = 175.3 . 
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Box 5a. N — 5, R = 0.23: the function 

j(t) := (2,r)3/2|„5^_^^^^(t)| Qnc has 7(0) = 

22.27..., 7(0.5) = S.ieO..., 7(1) = 2.997..., 
7(1.5) = 1.102..., 7(2) = 0.4055..., 7(4) = 
7.428... X 10"^, 7(8) = 2.491... X 10"^, 
7(10) = 4.564... X 10 
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Box 5b. JV = 5, fi = 0.23: the function 
'D3(t) = ||u^(t)||3. One has 153(0) = 154.3..., 
-03(0.5) = 56.97..., 1)3(1) = 20.90..., 
-03(1.5) = 7.646..., 153(2) = 2.810..., 
•03(4) = 0.05146..., -03(8) = 1.726... X lO"'"^, 
1^3(10) = 3.162... X 10"'' . 
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Box 5c. AT — 5, R — 0.23: the function 

£3(4). One has £3(0) = 0, £3(0.33) = 
1.613... X 10"^, £3(2) = 2.061... X lO"'', 
£3(3) = 8.234 X 10"^", £3(4) = 1.182... X 
10-11 . 
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Box 5d. N — 5, R = 0.23: the function 
3^3(t). This appears to be globally defined, 
and vanishing at +00. One has 3^3(0) — 
0, 3^3(1) = 1.004..., 3^3(1.5) = 2.609..., 
313(2) = 3.014..., K3(4) = 0.7907..., 3^3(8) = 
0.01580..., 3^3(10) = 2.141... X 10"^ . 
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Box 6a. N — 5, R — 0.24: the function 


7(t) := (2,r)3/2[„5^^^_^j(t)|. One has 7(0) = 


22.27..., 7(0.5) = 8.157..., 7(1) = 2.996..., 


7(1.5) = 1.101..., 7(2) = 0.4053..., 7(3) = 


0.05485..., 7(4) = 7.424... x 10"^ . 
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Box 6b. JV = 5. fl = 0.24: the function 


-03(4) = ||u'"(t)||3. One has I>3(0) = 154.3..., 


-03(0.5) = 58.08..., -03(1) = 20.90..., 


153(1.5) = 7.643..., 153(2) = 2.808..., 


153(3) = 0.3800, 153(4) = 0.05143... . 
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Box 6c. iV — 5, i? — 0.24: the function 


£3(t). One has £3(0) = 0, £3(0.33) = 2.0827X 


10"^..., £3(1) = 1.043 X 10~*..., £3(2) = 


2.664... X 10~^, £3(4) = 1.591... X 10~11 . 
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Box 6d. N = 5. R = 0.24: the function 


X.3(t). This diverges for t ^ T^ = 3.332.... 


One has 3^3(0) = 0, 3^3(0.5) = 0.06348..., 


3?3(1) = 2.126..., ^3(2) = 10.89..., 3^3(3) = 


33.31..., 3^3(3.33) = 4520.7... . 
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Box 7a. JV = 5, -B = 0.12: tho function 
7(t) := (2^)3/2|„l^_^_^j{t)|. One has 7(0) = 

22.27..., 7(0.5) = 8.i84..., 7(1) = 3.009..., 
7(1.5) = 1.107..., 7(2) = 0.4072..., 7(4) = 
7.459... X 10"^, 7(8) = 2.502... x 10~®, 
7(10) = 4.583... X 10"^ . 




12 3 4 

N = 5, R = 0.12: the function 

||ti^(t)||3. One has I>3 (0) = 154.3..., 

= 57.10..., Xi3(l) = 20.88..., 

= 7.672..., 113(2) = 2.821..., 

0.05168..., 153(8) = 1.733... X 10"^, 



Box 7b 

Vsit) : = 

©3(0.5) 

I'3(l-5) 

1^3(4) 

^3(10) = 3.175... X 10 
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Box 7c. 


iV = 5, fl = 0.12: tho function 


£3{t). One 


has £3(0) = 0, £3(0.33) = 


3.259... X 10" 


"^, £3(1) = 1.634... X 10^'^, 


£3(2) = 4.126 


... X 10""^, £3(4) = 1.266... X 


10-" . 
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Box 7d. JV = 5, ii = 0.12: the function 
3i3(t). This appears to be globally defined, 
and vanishing at +00. One has 3^3(0) — 0, 
313(1) = 2.237... X 10""'', 313(1.4) = 2.784... x 
10"*, 3^3(2) = 2.171... X IQ-*, 5:3(4) = 
3.410... X 10-^, 3^3(10) = 8.431... X 10"** . 
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Box 8a. Galerkin method with G as 

in 15.291 . for R = 0.12: the function 
7(t) := (27vf/^\uf^_^_g^(t)\. One has 7(0) = 

22.27..., 7(0.5) = 8.184..., 7(1) = 3.009..., 
7(1.5) = 1.107..., 7(2) = 0.4072..., 7(4) = 
7.459... X 10"^, 7(8) = 2.506... X 10"'', 
7(10) = 4.562... X 10""** . 
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Box 8b. Galerkin method tor R = 0.12: 



the function 'V>^(t) 
has 63(0) = 154.3.. 
©3(1) = 20.88..., 
153(2) = 2.821..., 
53(8) = 1.736... X 10" 
lO-'' . 



:= ||u''(t)||3. One 
, 63(0.6) = 57.10..., 
•63(1.5) = 7.672..., 
63(4) = 0.05168..., 
■^, 63(10) = 3.811... X 
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Box 8c. Galerkin method for R = 0.12: 
the function e^(t). One has £^3(0) — 0, 
^3(0. 18) = 1.431..., ?3(1) = 0.01121...^ 
^3(2) = 4.492... X 10"'', ?3(4) = 5.590... X 
10-^ . 




2 4 6 8 10 

Box 8d. _Galerkin method for R — 0.12: the 
function 3^3 (t). This appears to be globally 
defined, and vanishing at +00. One has 
S3(0) =_0, ^3(1) = 2.964..., ^3(1. 5) = 
4.022..., 3^3(3) = 1.6J9..., »3(4) = 0.6665..., 
3:3(6) = 0.09327..., 3^3(8) = 0.01267, $3(10) 
= 1.716... X 10"^. 
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